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ABSTRACT 

Context. Before the launch of the Swift satellite, the majority of the gamma-ray burst (GRB) afterglows for which Lya was redshifted 
into the observable spectrum showed evidence for a damped Lyo' absorber. This small sample indicated that GRBs explode either in 
galaxies, or regions within them, having high neutral hydrogen column densities. 

Aims. To increase the spectroscopic sample of GRBs with z> 2 and hence establish the N(ll l) distribution along GRB lines-of-sight. 
Methods. We have obtained six z > 2 GRB afterglow spectra and fitted the Lyo- absorption line in each case to determine ^(H l). This 
has been complemented with 12 other Swift Nifi l) values from the literature. 

Results. We show that the peak of the GRB A'(Hl) distribution is qualitatively consistent with a model where GRBs originate in 
Galactic-like molecular clouds. However, a systematic difference, in particular an excess of low column-density systems compared 
to the predictions, indicates that selection effects and conditions within the cloud (e.g. strong ionization) influence the observed 
A'(Hl) range. We also report the discovery of Lyce emission from the GRB 060714 host, corresponding to a star-formation rate of 
approximately 0.8 Mq yr"' . Finally, we present accurate redshifts of the six bursts: z = 3.240±0.001 (GRB 0503 19), z = 2. 198 ± 0.002 
(GRB 050922C), z = 3.221 + 0.001 (GRB 060526), z = 3.425 ± 0.002 (GRB 060707), z = 2.711 ± 0.001 (GRB 060714) and z = 
3.686 + 0.002 (GRB 060906). 

Key words, gamma rays: bursts - galaxies: high-redshift - galaxies: abundances - dust, extinction 



1. Introduction 

In just 18 months, the Swift satellite (Gehrels et al. 120041 has 
already contributed considerably to the progress of gamma-ray 
burst (GRB) science. It detects roughly two GRBs/week and 
transmits their accurate localizations (error radius frequently 
less than 5") to the ground within minutes for rapid follow-up 
observations. The GRB redshift distribution found by Swift is 
very different from that of the pie-Swift sample, being skewed 
to much higher redshifts (Berger et al. 120051 Jakobsson et 
al. |f2006a 2006b Daigne et al. 20Q6| Le & Dermer 2006). 
Approximately 70% of the Swift bursts are found to be located 
at z > 2, while the corresponding fraction was only 20% for 
pre-Swift bursts. 

It is thus much more common for the Lya line in Swift bursts 
to be redshifted redward of the atmospheric cutoff in optical 
spectra. This provides us with the opportunity to investigate the 
GRB host neutral hydrogen column density distribution. In par- 
ticular we can compare it to the damped Lya absorbers (DLAs) 
seen in absorption against QSO spectra (see Wolfe et al. 2005 for 
a recent review), defined as displaying A^(Hi) > 2 x 10^° cm"^ 



(logA^(Hi) > 20.3). Most of the neutral gas in the Universe in 
the redshift interval < z < 5 is in DLAs, providing the fuel for 
star formation at these epochs. Given that long-duration GRBs 
are known to have massive stellar progenitors (e.g. Hjorth et al. 
2003b, Malesani et al. 2004 Fruchter et al. 2006 1, GRB-DLAs 
provide valuable information on the sites of active star formation 
in the high-redshift Universe. In the pie-Swift era, only seven 
GRB A^(H i) column density measurements were obtained, with 
six being classified as DLAs (Vreeswijk et al. 120041 and refer- 
ences therein). 

In this Letter we present optical spectroscopy of six GRBs, 
focusing on the measurement of the A^(H i) column density. We 
also report the detection of Lya emission from one of the bursts 
(GRB 0607 14). We then discuss how the observed A^(Hi) col- 
umn density distribution from Swift bursts can be understood in 
terms of the GRB environment. 

2. Observations 

GRBs 050319, 050922C, 060526, 060707, 060714 and 060906 
are all long-duration bursts detected by the Swift satellite. Each 
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Table 1. A log of the follow-up spectroscopic observations for 
the six bursts presented in the paper. Af is the time from the onset 
of the burst. 



GRB 


Tel/instrument/grism 


Exposure 
time [s] 


A? 
[days] 


Spectral 
res. [A] 


050319 


NOT/ALFOSC/#4 


3 X 2400 


1.50 


7 


050922C 


NOT/ALFOSC/#4 


2400 


0.05 


5 


060526 


VLT/FORS 1/600 V 


900-1-1800 


0.44 


5 


060707 


VLT/FORS 1/300 V 


3x 1800 


1.44 


10 


060714 


VLT/FORS 1/300 V 


3x 1800 


0.50 


8 


060906 


VLT/FORS 1/300 V 


600 


0.05 


8 



burst was localized by the Burst Alert Telescope, each posi- 
tion refined by the X-Ray Telescope and a subsequent opti- 
cal afterglow (OA) was detected in every case (Rykoff" et al. 
l2005all2U05F ; Campana et al. '2006a' de Ugarte Postigo l2006l 
Krimm et al. 12006; Cenko et al. 2006b). 

Using the Nordic Optical Telescope (NOT), we obtained 
specti -a of the OA of GRBs 05 0319 and 050922C (Fynbo et al. 
EU05l Jakobsson et al. l2005b> . The data were acquired with the 
ALFOSC instrument with a r/3 wide slit and a grism with a 
wavelength coverage from 3700-9 100 A. The FORSl instru- 
ment on the Very Large Telescope (VLT) was used to obtain OA 
spectra of GRBs 060526, 060707, 060714 and 060906 (Thone 
et al. in prep; Jakobsson et al. 2006c 2006d Vreeswijk et al. 
12006 1. A r/0 wide slit was used and grisms with a wavelength 
coverage from 4700-7100A (GRB 060526) and 3600-9000 A. 
The details of our observations are given in Tabled] 

3. Results 

The spectrum of each OA displays a strong absorption line with 
usually clear damped wings. Blueward of this line, the flux 
drops substantially and exhibits the signature of the Lya forest. 
Associating the line with Lya, the rest of each spectrum was 
searched for additional features (>3cr) at the corresponding ap- 
proximate redshift. The results are presented in Table|2 

There is additional strong absorption in GRB050922C, 
corresponding to Lya at z 2.07. At this redshift we iden- 
tify four other features (Oi 1302/Sin 1304, Cii 1334/Cn* 1335, 
Sin 1526 and Civ 1548,1550) with an average redshift of z = 
2.075 ± 0.003. 

In Fig n we plot the normalized OA spectral region around 
Lya for each burst. Overplotted is a fit to the Lya absorption line 
yielding values reported in Table |2 Apart from GRB 060526, 
these values are well above the DLA definition. The redshifts 
deduced from the metal lines were found to be perfectly consis- 
tent with the DLA fine fits for GRBs 050922C and 060714. For 
the other four spectra, we note that slightly lower redshifts were 
used in order to provide the best fit to both the core and the red 
damped wing of the Lya profiles. However, adopting the metal 
line redshifts would make the A^(H i) only slightly smaller and 
within the reported Icr uncertainties. 

There is clearly a Lya emission in the centre of the 
GRB 060714 trough (Fig|2 see also Fig^, with an approxi- 
mate flux of 1.3 X 10"'^ ergs"' cm"^. In our assumed cosmol- 
ogy, (Qm, f^A, h) = (0.3, 0.7, 0.7), this coiTesponds to a lumi- 
nosity of 7.9 X lO"*' erg s"'. We can use this result to derive the 
star-formation rate (SFR), assuming that a Lya luminosity of 
10"*^ erg s"' c orresp onds to a SFR of 1 Mq yr"' (Kennicutt 199^ 
Cowie & Hu [T^I I. The Lya SFR in the GRB 060714 host is 
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Fig. 1. The normalized OA spectra centred on the Lya absorp- 
tion lines at the GRB host galaxy redshifts. A neutral hydrogen 
column density fit to the damped Lya line is shown with a solid 
line in each panel, while the Icr eiTors are shown with dashed 
lines. 
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Table 2. The redshift, neutral hydrogen column density and absorption lines identified for each afterglow spectrum. The number in 
square brackets after the redshift indicates the number of lines used to calculate it. 



GRB 




logMH I) 


050319 
050922C 


3.240 + 0.001 [3] 
2.198 + 0.002 [8] 


20.90 ± 
21.55 ± 


0.20 
0.10 


060526 


3.221 ± 0.001 [12] 


20.00 ± 


0.15 


060707 


3.425 + 0.002 [8] 


21.00 ± 


0.20 


060714 


2.711 +0.001 [36] 


21.80 ± 


0.10 


060906 


3.686 ± 0.002 [5] 


21.85 ± 


0.10 



Absorption features 



Lyff, Sill 1260, Oi 1302, Cii 1334/Cii* 1335, Civ 1548,1550 

Lyff, Sin 1260, Sin* 1264, Oi 1302/Siii 1304, Cii 1334/Cn* 1335, Siiv 1393,1402, 

Sin 1526, Sin* 1533, Civ 1548,1550, Feii 1608, Alii 1670 

Sin 1190,1193, Sim 1206, Lya, Sin 1260, Oi 1302, Sin 1304, Cn 1334, 

Cii* 1335, Sin 1526, Sin* 1533, Fen 1608, Aln 1670 

Ly/J, Cn 1036/Cn* 1037, Nn 1083, Sim 1206, Lyo-, Sin 1260, Oi 1302/Sin 1304, 
Cn 1334/Cii* 1335, Sin 1526, Civ 1548,1550, Aln 1670 

Lyy, Cm 977, Sin 989, Ly/3, Ovi 1031,1037, Cn 1036, Cn* 1037, Nn 1083, Fen 1144, 
Sin 1190,1193, Lya, Sn 1250,1253, Sin 1260, Sin* 1264, Oi 1302/Sin 1304, Siii* 1309, 
Cn 1334/Cii* 1335, Nin 1317,1370, Siiv 1393,1402, Sin 1526, Sin* 1533, Civ 1548,1550, 
Fen 1608, Aln 1670, Nin 1709,1741, Sin 1808, Aim 1854,1862, Znn 2026,2062, 
Fen 2260,2344,2374,2382, Fen** 2328,2349, Fen*** 2359, Fen* 2365 
Ly/J, Cn 1036/Cii* 1037, Lya, Sin 1260, Oi 1302/Sin 1304, 
Cn 1334/Cii* 1335, Siiv 1393,1402, Sin 1526 




Fig. 2. The two-dimensional GRB 060714 OA spectrum centred 
on the Lya absorption line. Lya in emission is clearly extended 
and centred in the trough. 



thus ~0.8Moyr which is within the range of values found 
for o ther GR B hosts (e.g. Fynbo et al.|2P02 2003 Jakobsson 
et al. l200 5aV We note that this value has not been corrected for 
host extinction, and is therefore a strict lower limit to the actual 
SFR. Fynbo et al. l2003!l find that Lya emission is much more 
frequent among pre-5w//f GRB host galaxies than among the 
Lyman-break galaxies at similar redshifts and that a low metal- 
licity preference for GRBs could be the explanation (see also 
Vink & de Koter 2005 , Fruchter et al. 2006 , Priddey et al. 2006 ). 
We are currently undertaking a survey of GRB host galaxies that 
will constrain the fraction of Lya emitters among Swift GRB 
hosts. 

Although there are many metal lines observed in the GRB 
spectra, most of them are saturated. However, we have identified 
a couple of lines, corresponding to elements that are known to 
deplete negligibly on dust. They fall well outside the Lya forest 
and are thus probably unblended. In addition, these transition 
lines are very weak, implying they are the least saturated ones: 
Sin 1808 and Znii 2026 in GRB 060714. In the optically thin 
limit approximation, their equivalent widths correspond to the 
following metalHcity limits: [Si/H] > -1.35 and [Zn/H] > -1.0. 



4. Discussion 

In Table |3] we have compiled all available H i column density 
measurements for z > 2 Swift bursts (currently displaying a me- 
dian redshift of 3.4). The bottom panel of Fig.|3]shows a com- 
parison between these H i column densities and those of the pre- 
Swift sample (Vreeswijk et al. i2004j . The former currently has 



Table 3. A list of all Swift GRBs with z > 2 known to date 
(1 October 2006). A^(Hi) is derived from optical spectroscopy. 
References are given in order for the redshift and the H i col- 
umn density (information is not currently available for bursts 
marked with — ). (1) This work; (2) Watson et al. ( 2006a) ; (3) 
Be rger et al. (200^; (4) Berger et al. (|2005); (5) Starling et 
al. i2005t ; (6) Ledoux et al. ( l2U05t ; (7) Ledoux et al. (in prep); 
(8) Kawaiet al. (|2006); (9) Foley et al. (2005 1; (10) Quimby 
et al. (|2M5^; (11) Piranomonte et al. (2006); (12) Prochaska 
et al. ( I2006bi ; (13 ) J. X. Prochaska (private communication); 
(14) Fynbo et al (I2006b> ; (15) Cucchiara et al. (|2 006ai; (16) 
Cucch iara et al. (|2006b); (17) Berger et al. f2006a); (18) Price 
(l2006b ; (19) Cento et al. (2006a); (20) Ferrero et al. (2006); (21) 
Ledoux et al. (2006|l; (22) Smette et al. (in prep); (23) Rol et al. 
FmO^; (24) D'Eliaetal. ( l2006i ; (25) Fugazza et al. (in prep); 
(26) Fynbo et al. Il2()()6a> . 



GRB 


7 


logA'(Hl) 


References 


050319 


3.24 


20.9 ± 0.2 


(1)(1) 


050401 


2.90 


22.6 ± 0.3 


(2) (2) 


050505 


4.27 


22.1 ±0.1 


(3) (3) 


050603 


2.82 




(4) 


050730 


3.97 


22.1 ±0.1 


(5) (5) 


050820A 


2.61 


21.1 ±0.1 


(6) (7) 


050904 


6.30 


21.3 


(8) (8) 


050908 


3.34 


19.2 


(9) (9) 


050922C 


2.20 


21.6 ±0.1 


(1)(1) 


051109 


2.35 




(10) 


060115 


3.53 




(11) 


060124 


2.30 


<20.3 


(12) (13) 


060206 


4.05 


20.9 ±0.1 


(14) (14) 


060210 


3.91 


21.7 ±0.2 


(15) (16) 


060223 


4.41 




(17) 


06051 OB 


4.94 




(18) 


060522 


5.11 


20.5 ± 0.5 


(19) (19) 


060526 


3.22 


20.0 ± 0.2 


(1)(1) 


060605 


3.71 




(20) 


060607 


3.08 


<19.5 


(21) (22) 


060707 


3.43 


21.0 ±0.2 


(1)(1) 


060714 


2.71 


21.8 ±0.1 


(1)(1) 


060906 


3.69 


21.9±0.1 


(1)(1) 


060908 


2.43 




(23) 


060926 


3.21 


22.7 ±0.1 


(24) (25) 


060927 


5.6 




(26) 
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a median logA^(Hi) of 21.6, a bit higher than the latter (21.3). 
However, the small size of the pre-5w//f sample does not allow 
us to reject, with any degree of confidence, the null hypothesis 
that the two samples are drawn from the same distribution, i.e. a 
two-sample Kolmogorov-Smirnov (KS) test results in a 90% sig- 
nificance. The fraction of DLAs is also similar in both samples 
(~80%), supporting the conclusion that GRB absorption systems 
show exceptionally high column densities of gas when com- 
pared to DLA systems observed in the lines-of-sight to QSOs 
(Vreeswijk et al. .2004rf . Indeed, a two-sample KS test indicates 
there is a less than 10"' probability that GRB-DLAs and QSO- 
DLAs (Prochaska et al. 2005 1 are drawn from the same popula- 
tion. This result presumably reflects the fact that GRBs occur in 
star-forming regions within their hosts, whilst QSOs select more 
random lines-of-sight through intervening galaxies. 

In the bottom panel of Fig.l^we also compare the observed 
GRB A^(H i) distribution to the expected column density distri- 
bution for bursts in Galactic-like molecular clouds (Reichart & 
Price UJW\ hereafter RP02). The model (solid, non-filled his- 
togram) is mass-weighted and corrected for geometrical effects, 
i.e. the clouds being centrally condensed, and the GRB loca- 
tion within a cloud (not behind it). Qualitatively, especially in 
terms of the location of the peak of the distribution, the match 
is relatively good. However, a KS test shows that the model and 
observed A^(H i) distributions are inconsistent with being drawn 
from the same population at the 3cr level, with a clear overabun- 
dance of low-A^(H i) detections compared to the model. A possi- 
bility is that GRB host galaxies tend to have lower A^(H i) clouds 
than the Milky Way (MW) and/or GRBs are more likely to oc- 
cur in lower A^(H i) clouds. Some evidence for the former exists 
(Rosolowsky 2005 1, with the mass distribution of clouds in the 
LMC (which is more similar to GRB hosts than the MW) having 
a marginally steeper distribution than the inner disk of the MW. 

A few other eff'ects can give rise to the excessive detections 
of logA^(Hi) < 21.0. The GRB progenitor and nearby mas- 
sive stars may (partially) ionize their local environments; there 
is clearly a trend for GRBs with the smallest A^(H i) to exhibit 
very weak low-ionization lines (e.g. Si ii and C n) while display- 
ing sti-ong absorption of e.g. Si iv and C iv (GRB 021004: M0ller 
et al. 2002 GRB 050908: Prochaska et al. 2006a GRB 060124: 
Prochaska et al. 2 006bj . This may support the high-ionization 
scenario for the low-A^(H i) systems, although implicit is the as- 
sumption that low- and high-ionization species trace the same 
phase. An alternative scenario is that some GRBs are formed by 
massive runaway stars (e.g. Hammer et al. 2006 i; these would 
explode in regions of relatively low A^(H i). Finally, it has been 
su ggeste d (e.g. White et al. 1999, Sugitani et al. 120021 Miao et 
al. '2006 1 that a significant part of star formation in molecular 
clouds takes place at their outer edges. The triggering of this star 
formation also blows the cloud open, resulting in these stars hav- 
ing (largely) unimpeded view outside the cloud and presumably 
a low A^(H i) column. 

The drop-off in detections above log A^(H i) > 22.0 could be 
due to selection effects, i.e. the sample might be incomplete due 
to high extinction. GRBs most likely do not burn through all of 
the dust in their molecular clouds (e.g. fig. 4 in RP02), so some 
high-A^(H i) systems are likely dimmed by dust and consequently 
missed in the sample in Table|3(see also Savaglio et al. 2003 
Vergani et al. 2004 1. In the top panel of Fig. we have plotted 
the corresponding afterglow 7?-band magnitudes (UV restframe) 
interpolated to a restframe epoch of 12 hr and redshifted to z = 3. 
The majority of the bursts have a well-sampled light curve, al- 
lowing us to determine the magnitudes via spline interpolation 
(see appendix A in Kann et al. .2006 for a detailed description 




19 20 21 22 25 

log N(HI) 



Fig. 3. Bottom: Histogram of the H i column density measured 
in GRBs for which the redshift was large enough to detect Lya. 
Bursts detected by Swift are filled while the ^r&Swift sample 
is hashed (Vreeswijk et al. I2004> . The current range of known 
GRB Hi systems is 19.2 < logA^(Hi) < 22.7. For comparison, 
the QSO-DLAs (Prochaska et al. I2005t are overplotted as the 
dashed histogram (normalized with respect to the GRB obser- 
vations). The solid, non-filled histogram shows the GRB A^(H i) 
prediction by RP02. Top: The corresponding afterglow 7?-band 
magnitudes, corrected for Galactic extinction (Schlegel et al. 
1998), interpolated to a common epoch of 12hr (in the source 
restframe) and redshifted to z = 3. GRBs 060210 and 060522 are 
omitted due to the lack of information on their spectral energy 
distribution; hence corrections could not be made to account for 
Lya absorption. Where the uncertainty on the H i column is not 
reported in the literature, a value of 0.2 has been plotted. 

of the method used). For the most recent bursts with a relatively 
poor data sampling, we have applied extrapolation and adopted 
large conservative error bars. There is only a tentative evidence 
of an anti-correlation between the OA flux and higher A^(H i) 
at the 2cr level as derived from Spearman's rank correlation. It 
should be noted that intrinsic GRB /?-band afterglow light curves 
show evidence for a bimodal luminosity distribution (Kann et al. 
2006; Liang & Zhang 2006; Nardini et al. 2006 1, that increases 
the scatter in the top panel of Fig.|3 Excluding the upper limits, 
there is indeed an indication of bimodality in the data. 

The characteristics of the dust-to-gas ratio in the clouds are 
discernible from the range of more than three orders of magni- 
tude in the H i column densities (lower panel of Fig.|3}, suggest- 
ing a similar range in the dust column (e.g. Predehl & Schmitt 
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I1995> . If dust were entirely responsible for all of the observed 
variation in the /?-band (about five magnitudes; upper panel of 
Fig-E}, the range would, at most, be of that order, limiting it 
to a maximum of As a; 5 mag. At the Galactic dust-to-gas ratio, 
the observed As (UV restframe) corresponding to the maximum 
A^(Hi) value is around 80 magnitudes (e.g. Pei ll992l Diplas & 
Savage I994l, implying that the As/A^(Hi) ratio in these ab- 
sorbers must be at most 7% of the Galactic value (see also Hjorth 
et al. 2003a). If some of the variation in the flux is due to other 
factors, the ratio drops even further. 

Given that it is thought that most of the star formation in the 
Universe occurs in molecular clouds, it seems logical that long- 
duration GRBs originate in such regions (Fig.|3}. However, there 
is a systematic difference between the model and observations 
that presumably can be clarified as a result of a combination of 
selection effects (dust extinction) and intense ionization in the 
GRB environment. Combining the A^(Hi) values derived from 
optical spectroscopy with the metal column densities from soft 
X-ray absorption would conceivably shed further light on these 
issues. A detailed analysis of this is presented in a subsequent 
paper (Watson et al. 2006b i. 
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